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ABSTRACT 

We present results from a spectral analysis of a sample of high-redshift (z > 3) X- 
ray selected AGN in the 4 Ms Chandra Deep Field South (CDF-S), the deepest X- 
ray survey to date. The sample is selected using the most recent spectroscopic and 
photometric information available in this field. It consists of 34 sources with median 
redshift z — 3.7, 80 median net counts in the 0.5-7 keV band and median rest-frame 
absorption-corrected luminosity L 2 - io kcv ~ 1.5x 10 44 ergs _1 . Spectral analysis for the 
full sample is presented and the intrinsic column density distribution, corrected for 
observational biases using spectral simulations, is compared with the expectations of 
X-ray background (XRB) synthesis models. We find that s» 57 per cent of the sources 
are highly obscured (Njj > 10 23 cm~ 2 ). Source number counts in the 0.5 — 2keV band 
down to flux F .5_2koV ~ 4 x 10 _17 erg s _1 cm~ 2 are also presented. Our results are 
consistent with a decline of the AGN space density at z > 3 and suggest that, at those 
redshifts, the AGN obscured fraction is in agreement with the expectations of XRB 
synthesis models. 
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1 INTRODUCTION 

The discovery of Super Massive Black Holes ( SMBH) at z > 
6 with masses of the order of 10 9 M 
20111. the redshift record z 



7.085; 



Mortlock et all 
, IWillott et all |2003t 
Fan et al. l l2006al lb HWillott et al. 2009) implies that SMBH 
were already in place when galaxies were only a few hun- 



dred Myr old. The formation of these massive objects 
and the physical conditions allowing a rapid and early 
growth are challenging problems. The origin of the SMBH 
seeds is also uncertain. Two classes of formation mod- 
els are currently under discussion: light seeds, perhaps 
the r e mnants of Pop HI sta r s (~ 1 2 M^; iMadau fc Reed 
20011; lAlvarez Wise fc Abell 120091: [Volonteri fc Begelmanl 
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201C . see also Whalen fc Frvei 12011 ). possibly evolved as 
a "quasi-star" l|Begelman |2010| ). and heavy seeds, formed by 
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monolithic colla pse of massive gas clouds (10 4_5 M o ; e.g. 
iBall et alJlioilh . Whatever the nature of the seed is, contin- 
uous Eddington-limited or even super-Eddington accretion 
onto these objects is required to reach the observed masses. 

Such a high accretion rate requires a large amount of 
accreting material which, at the same tim e, may act as 
a scr een against the emitted radiation (e.g. iHopkins et al.l 



2008). Supporting observational evidence was found for 
an increasing fraction of obscured Active Galactic Nuclei 
(AGN) from z = to z sa 2 - 3 dLa Franca et al.l [2005 : 



iBallantvne et al.ll2006l;lTreister fc UrrvteOOfllHasingerlhoOc 
but see also iGilli et al.l l20ld . who dispute these results) 
At these redshifts, a large amount of gas is available in 
the host galaxy. The nuclear activity could be triggered 
by the gas falling towards the centre, when a "wet" (i.e. 
gas-rich) merging event occurs. After an initial obscured 
(Nh > 10 22 cm -2 ) accretion phase, which in some models 
can be Compton-thick (Nh > 10 24 cm~ 2 ), the emitted radi- 
ation sweeps out the gas, revealing itself and acting as a neg- 
ative feedb ack towards the star formation in t he host galaxy 



(e.g.lBrusa et al. 2005; Alexander et al 



2010i; ISturm et aljfe oil: Mai olino et al 



201fj; lFeruglio et al.l 
2012T ). The evolution 



of obscuration, in particular at high redshift, when most of 
the mass-accretion was taking place, is a fundamental ob- 
servable to probe this framework. 

The AGN comoving space density in the optical and 
X-ray band is not constant over time: there are many more 
AGN at z > 1 than in the local Universe (e.g. lHasinger et al.l 
2005). The evolution of the AGN population can be de- 
scribed by a phenomenological model known as Luminos- 
ity Dependent Density Evolution model (LDDE model; 
e.g. lHasinger. Mivaii fc Schmidtl2005l ; lBrandt fc Alexander! 
2010): according to this model, the comoving space den- 
sity of more luminous AGN (Lx > 10 45 ergs _1 ) peaks at 
z ~ 2 — 3 and then decreases exponentially to z ~ 6; for less 
luminous ones (L x < 10 45 ergs _1 ) the peak arises at more 
recent times (z ~ 1 — 1.5) and then their density slowly 
decreases till the highest redshifts so far probed (z ~ 3). 
Therefore, more massive black holes formed and grew ear- 
lier than smaller ones. This sort of "cos mic downsizing" is 
similar to that found for galaxies (e.g . ICowie et alj 1 19991 : 
iThomas et aHl2005l ; [Pamen et al.ll2009l ). 

This scenario is, however, far from being complete. Par- 
ticularly, the behaviour at higher redshifts (z ^ 3) is not 
constrained (especially at Lx < 10 44 ergs _1 ), implying se- 
vere uncertainties in AGN evolution during the early ages 
of the Universe. Furthermore, the space density and the role 
of heavily obscured AGN are largely unknown beyond the 
local Universe. In order to place constraints on the obscured 
AGN fraction and its evolution history, observables, such as 
the spectral properties and the number counts, are required 
to allow comparisons with evolutionary theoretical models. 

Optical surveys are not suitable for these tasks since 
they are severely biased against obscuration. Deep X-ray 
surveys are required, because they are less affected by ob- 
scuration and can reach fluxes faint enough to sample the 
bulk of the high-redshift AGN population. Tricha s et al.l 
|2012l ) presented the largest (78 objects) sample of z > 3 
X-ray selected AGN, using data from the Chandra Multi- 
wavelength Project (ChaMP) survey. A few studies have 
been performed on high-redshift AGN exploiting the deep 
surveys in the Cosmological Evolution Survey (COSMOS) 



field carried o ut with XMM-New ton (|Brusa et al.l 120091 ^ 
and Chandra (|Civano et a.1.1 l201ll ) , and they were limited 
to relatively bright X-ray fluxes (Jo.5-2kcV 

> 10~ 15 and 

3 x 10 -16 ergcm _2 s _1 , respectively) and high 2 — lOkeV 
luminosities (Z/2-iokcV > 10 44 ' 2 ergs _1 and 10 43 ' 5 ergs -1 , 
respectively). Only recently, t hanks to the 4 Ms Chan- 
dra Deep Field South CCDF-S. Ixiie et alJlioTH ). the deep- 
est X-ray survey to date, very faint fluxes (down to F ~ 
3.2 x 10~ 17 /9.1 x 10" 18 /5.5 x 10~ 17 ergem - ^ -1 in the 
0.5 - 8/0.5 - 2/2 - 8keV band) could be reached. Early 
results on z > 3 p opulation from the 4 Ms CDF-S da t a have 
been presented in iFiore et al] (|2012T l. iLehmer et ail (|2012T ) 
also presented the source counts of the z > 3 population 
as part of a more comprehensive work on logN-logS of the 
X-ray population. 

In this work, we analyse data from the 4 Ms CDF-S. In 
particular, we have selected a sample of "bona fide" high- 
redshift (z > 3) AGN using the most recent spectroscopic 
and photometric information available, and performed an 
X-ray spectral study to understand their properties and, 
possibly, constrain the early phases of AGN evolution. The 
spectral properties, column-density distribution and number 
counts of the sample are derived, discussed and compared 
with previous works and model predictions. 

Throughout this work we assume a Sl m = 0-3, S1a = 0. 7 
and H = 70 km s^Mpc" 1 Universe (ISpergeletal .1120031 ). 
The errors on the spectral parameters correspond to t he 90 
per c ent confidence level for one parameter of interest |Avnil 

Eizi). 



2 THE z > 3 AGN SAMPLE 



We assembled a sample of X-ray 
z > 3 in the 4 Ms CDF-S 
both spectroscopic and photometric 



scopic redshi fts 



for 



selected AGN at 
on the basis of 
redshifts. Spectro- 
> 3 AGN were collected from 



ISzokqlv et~aH (120041). IVanzella et~aT1 j|2008h. iPopesso et~ail 
<|2009|). ISilyerman et ail (|2010l ), IVanzella et alj (|2010T ). 
IWuvts et al.l |2009l ) and E. Vanzella (private commu- 
nicatiorQ) . A qua lity flag ( " secure /insecure") is associ- 
ated to them in IXue et ail (1201 ll). excep t for spectro- 
scopic redshifts from Vanzella et alj l|2010l ), IWuvts et al.l 
(2009) and E. Vanzella (p ri vate communication) , which 
are not used by IXue et ail (1201 ll ) and for which we re- 
trieved the original quality information . Phot ometric red- 
shifts were gathered from iLuo et al.1 (|2010|. which also 
provides seco n dary so lutions), 



. .Cardamone et al.1 d2O10Tl . 
iRafferty et ail (I2OIII), ISantini et al.1 (120091) IWuvts et al.1 



Kattcrty • 

[ 20081) . iTavlor et al.1 (|2009l ) and lDahlen et al.1 (|2010l ). 

We first cross- matched the X-ray positions of all sources 
in th e 4 Ms CDF-S main catalogue (740 sources; IXue et ail 
2011) with the positions of the sour c es in the catalogues 
which are not collected by IXue et al.1 (|201lD using a 1 arc- 
sec matching radius. Considering the ot her catalogues, w e 
assumed the counterparts reported by IXue et alT l|201 ll ). 



1 These spectroscopic redshifts are based on a re-analysis of Vis- 
ible MultiObject Spectrograph (VIMOS) and FOcal Reducer and 
low dispersion Spectrograph (FORS2) spectra. 
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Figure 1. Selection criteria applied to the 96 pre-selected sources (which have at least one spectroscopic or photometric redshift z > 3) 
to be included in the high-redshift sample, visualised as a flow chart. The starting block has bold sides. 



which used a likelihood-ratio matching technique. We pre- 
selected a sample of high-redshift AGN candidates, includ- 
ing all sources which have a spectroscopic or photometric 
redshift z > 3 in at least one of the considered catalogues (96 
sources). Given the large amount of available information, 
we built up a set of criteria to select the final high-redshift 
sample. A source is included in the sample if it fulfils one of 
the following requirements: 

(1) A spectroscopic redshift flagged as "secure" z speCta > 
3 is provided. 

(2) A spectroscopic redshift flagged as "insecure" 
z sp e C ,i > 3 is available and it is upgraded to "secure" by 
re-analysing the optical spectrum, checking the quality and 
reliability of the identified spectral features. 

(3) N phot , >3 - N pho t,<3 > 3, where N phott>3 is the num- 
ber of photometric redshifts larger than 3 and N p hot,<3 is 
the number of photometric redshifts lower than 3, when no 
secure spectroscopic redshift and more than two photomet- 
ric redshifts are available. 

(4) When no secure spectroscopic redshift and two pho- 
tometric redshifts are available, both photometric redshifts 
are larger than 3. 



(5) When no secure spectroscopic redshift and only one 
photometric redshift is available, its la lower bound (pro- 
vided by the original catalogue) is larger than 3. 

The selection procedure is shown in Fig. [1] as a flow chart. 
We will discuss in § [4] the effects of relaxing the selection 
criteria on the sample size. 

The final sample consists of 34 sources (see Tab. [1} . A 
spectroscopic redshift is adopted for 15 sources (12 and 3 are 
included in the high-redshift sample on the basis of the first 
and second criterion, respectively) and a photometric red- 
shift is adopted for 19 sources (7/8/4 sources are included 
in the sample on the basis of the third/fourth/fifth crite- 
rion, respectively ). XID 392 has no redshift information in 
Xue et ah! (12011 ) a nd a photometric redshift z = 6.22 from 
Santini et alT l|2009). but no error on the redshift is provided. 
Therefore, we conservatively excluded this source from the 
sample. The sample presented i n this work is a fact or ~ 2 
smaller than that discussed in iLehmer et all (|2012T ). This 
difference will be discussed in more details in § [4j 

When more than one redshift per s ource is available , 
we almost entirely followed the choice of IXue et al.l |201l| y 
Briefly, the adopted redshift in this work is, in order of pref- 
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erence: secure (or upgraded to secure) s p ectro scopic redshift, 
photometric redshift from Luo et al.l l|2010l ). photometric 
redshift from lCardamone et al.l ( 20101 ). XID 331 has no red- 
shift i nformation in the catalogues collected by IXue et all 
(2011), but the inclusion of the other above-mentioned cat- 
alogues allows that source to be inc luded in the s a mple . We 
assumed photometric redshift from iDahlen et ail (120101) for 



10 



No spectroscopic redshift is reported in Xue et al. ( 20 111 ) 



for XID 235, 458, 528 and 386; redshift for these sources were 
obtained by E. Vanzella ( XID 2 35 and 458, pri vate commu- 
nicati on), IVanzella et al. | (|2010l . XID 528) and lWuvts et al.l 
( 2009, XID 386). A high quality flag is associated to all these 
redshift in the original works; therefore, we consider them 
as "secure". 

The redshift distribution of the sample is shown in 
Fig. [2] The mean and median redshifts are 4.0 and 3.7, re- 
spectively. Most of the sources with a spectroscopic redshift 
(13 out of 15) are at z < 4. We note the presence of 3 
sources at z > 5, for which we investigated in more de- 
tail the redshift information. Two of them (XID 139 and 
197) have only one phot ometric redshift (z — 5.729 and 
6.071, respectively, from iLuo et all [ioiol , which also pro- 
vides a secondary solution of z = 4.385 for XID 197), while 
XID 485 has two diff erent photometric redshifts (z = 7.62 
from ILuo et ah! |2010 | . with secondary solution z = 3.309, 
and z = 4.42 from lSantini et al.l l2009). These objects have 
no optical counterparts in the Great Observatories Origins 
Deep Survey (GOODS; XID 139 and 485) or Galaxy Evolu- 
tion from Morphologies and SEDs (GEMS; XID 197, which 
is not in the GOODS field) observations, while they are sig- 
nificantly detected by the Infrared Array Camera (IRAC) in 
the near-IR. Thus, a high-redshift (z > 3) solution appears 
probable. 



3 DATA ANALYSIS 

3.1 Spectral extraction procedure 

Chandra data products for each of the 54 observations of the 
CDF-S are publicly available in the Chandra Data Archive^. 
We retrieved all data products (in particular, the evt2 files; 
i.e., the event list filtered on the good time intervals and sta- 
tus bits), but those related to the first two observ ations (for 
which the focal-plane temperature was —110°; ILuo et al] 
l2008l ; IXue et alj|201ll ). and used them for spectral analysis. 
The sources were located in each of the 52 CDF-S images 
according to the coordina tes provided by the main-source 
catalogue (|Xue et al.ll201ll ). 

We selected the extraction regi ons for the source an d 
background using SAOImage DS9l (|jove fc Mandell lioolh . 
Despite the Chandra PSF distortion with the off-axis angle, 
we always chose circular regions, provided that most of the 
source counts are included in the selected regions (probably 
this choice slightly lowers the signal-to- noise ratio). Small 
radii are required for faint sources in order to include a low 
number of background counts and therefore to increase the 
signal-to-noise ratio. Larger extraction radii can be used for 
bright sources. Also the position of the source in the field 



Figure 2. Redshift distribution of all 34 sources (solid line) and 
sources with spectroscopic redshift (shaded histogram) in the 
sample. 



of view (f.o.v.) must be taken into account, because of the 
degradation of the PSF with the off-axis angle. For each 
source, we defined the same radius for the extraction re- 
gion in all observations, consideri ng the mean off-a xis angle, 
provided by the main catalogue |Xue et al.ll20Tll) . and the 
relation between it and the encircled energy radius. 

The off-axis angles for the sources in the sample are 
in the range of 1.42 — 10.1 arcmin. We chose the extraction 
radius for each source to obtain an acceptable trade-off be- 
tween the number of counts and the signal-to-noise ratio. 
As a general rule, for sources at off-axis angles < 4 arcmin 
we defined the extraction radius corresponding to an encir- 
cled energy fraction (EEF) of 90 per cent at 1.49 keV. For 
sources at a large off-axis angle (> 7 arcmin), we consid- 
ered radii in a range corresponding to 50-90 per cent EEF 
at 1.49 keV. For sources at intermediate off-axis angles (4-7 
arcmin) we tried different extraction radii and chose those 
which provided the best trade-off between source and back- 
ground counts. In a few cases we were forced to use smaller 
radii to avoid the contamination by sources close to the one 
we were considering. As a result, the extraction radii for 
the sample are in the range of 1.5 — 5arcsec. Regarding the 
background, we chose extraction regions as close as possible 
to the source, preferably on the same chip, in areas where 
no other source was detected. In particular, we excluded all 
the regions around detected sources with radii a factor ^1.2 
larger than the radii corresponding to 90% EEF at 6.4 keV 
(which are a factor > 1.5 larger than the radii corresponding 
to 90% EEF at 1.49 keV), in order to avoid contamination. 

We used CIAO v4.3jjto extract source and background 
spectra and to create the response files for each observation. 
In particular, we used the specextract tool which contains 
some improvements with respect to similar scripts in pre- 



2 http://cxc.harvard.edu/cda/ 

3 http:/ /hea-www. harvard.edu/RD/ds9/ 



http:/ / cxc.harvard.edu/ciao/ 
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Figure 3. Full-band (0.5 — 7keV) net counts distribution of all 
34 sources (solid line) and sources with spectroscopic redshift 
(shaded histogram) in the sample. 
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Figure 4. Best-fitting Nh as a function of the hard-band intrinsic 
luminosity for all the sources in the sample. Filled circles represent 
constrained values, while upper limits are reported as downward- 
pointing arrows. 



vious CIAO versions (e.g., it applies the energy- dependent 
aperture correction to the ARF file, accounting for the frac- 
tion of the PSF enclosed by a region). The spectra and 
response files of each source obtained for each observation 
were then added using tools provided in t he FTOOLS pack - 
a g43(i- e - mathpha, addarf and addrmf: iBlackburnl ll99ot ). 
weighted by the exposure time of the individual observa- 
tions. 

The distribution of full-band (0.5 — 7keV) net (i.e. back- 
ground subtracted) counts (Fig. [3} is peaked near 100 counts 
(the mean and median values are 275 and 80, respectively). 
As expected, the bulk of the objects has a very low number 
of net counts: ~ 74 per cent of sources in the sample have 
less than 200 net counts. In order to maximise the signal- 
to-noise ratio (SNR), in most cases we fitted the spectra 
in a narrower energy range (0.5 — 5keV), given that in the 
5 — 7 keV band the background is high. All the spectra were 
fitted using the Cash statistic <|Cashlll979l ) to estimate the 
best-fitting parameters. 

3.2 Spectral analysis 

The spec t ral fit ting was performed with XSPEC vl2.5.1rQ 
|Arnaudl [l996). We assumed an absorbed power-law, 
parametrized by XSPEC models powerlaw and zwabs, as the 
starting spectral model. Galactic absorption was included 
assuming a column d ensity in the range (6 — 9) x 10 19 cm -2 
jKalberla et al.ll2005l ). 

The mean photon index is (r) = 1.51 with dispersion 
a — 0.28 for the 8 brightest sources (> 250 net counts) and 
(r) = 1.63 with a = 0.89 for the whole sample. The aver- 
age photon index is slightly flatter than what is typically 

5 http: / /heasarc. nasa.gov / ftools / ftools_menu.html 

6 http:/ /heasarc. nasa.gov/xanadu/xspec/ 



measured in AGN samples (r = 1.8; e.g. iTozzi et al.ll2006l : 
iTurner et alll 19971 ). but we emphasize that we are sampling 
the rest-frame 2 — 20 keV band, and the likely presence of a 
reflection hump could flatten the observed X-ray spectrum. 
However, the measured value is statistically consistent with 
F = 1.8. Reasonable accuracy was obtained for the best- 
fitting parameters only for the 8 brightest sources (with a 
mean error on V and Nh of 14 and 45 per cent, respec- 
tively). The spectral parameters were not well constrained 
for almost all the other sources (error on F of 117 per cent; 
only upper limits on iV^ for most objects) because of the 
poor spectral quality. 

The difficulties in simultaneously constraining V and 
Nh for sources with few counts can be explained consid- 
ering that these two spectral parameters have a significant 
degree of degeneracy. The problem is further complicated by 
the redshifted absorption cut-off, which moves outside the 
Chandra band-pass at z > 3 even for significant degrees of 
obscuration, introducing further uncertainties. 

To better constrain the column density, we fixed the 
photon index to V = 1.8 and repeated the spectral analysis 
on all the sources. This photon index is a widely used value 
for sources with low counting statistics, being considered the 
typical slope for AGN power-law em ission in X-rays (e.g. 
ITurner et allll997l : ITozzi et all 120061 ) . The best-fitting pa- 
rameters, fluxes and intrinsic luminosities derived using this 
spectral model are reported in Tab. [1] Similar results are 
obtained fitting simultaneously the source and background 
spectra. All the X-ray spectra of the sources, fitted using 
this model, are shown in Appendix B. The effects of assum- 
ing a flatter photon index (T = 1.6) will be discussed in § 13.31 
and §H 

An anti-correlation between the column density and 
the AGN intrinsic luminosity is reported i n several works 
(e.g. ISteffen et all l2003t lUeda et alj 120031 . lHasingerl |2004 
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Table 1. Main information and spectral parameters of the 34 high-redshift AGN. 



XID 

(1) 


RA[J2000] 
(2) 


DEC[J2000] 
(3) 


z 

(4) 


r 

(5) 


C 
(6) 


N H 
(7) 


CSTAT/DOF 

(8) 


^0.5-2 kcV 

(9) 


^2-10keV 
(10) 


-^2-10 koV 
(11) 


COUNTS 
(12) 


27 


52.96054 


-27.87706 


4.385 


P8 


5 


«££ 


71.7/113 


2.45E-16 


1.76E-15 


2.80E44 


77 


100 


53.01658 


-27.74489 


3.877 


P8 


3 


oo+27 
38_ 24 


69.0/92 


1.48E-16 


<7.55E-16 


8.94E43 


67 


107 


53.01975 


-27.66267 


3.808 


P9 


1 


< 71 


51.6/62 


4.80E-16 


<1.40E-15 


1.52E44 


36 


132 


53.03071 


-27.82836 


3.528 


P8 


3 




69.0/65 


1.13E-16 


<3.39E-16 


3.12E43 


57 


139 


53.03496 


-27.67975 


5.729 


P8 


5 


96lf 5 


92.8/101 


5.71E-16 


3.64E-15 


1.02E45 


111 


150 


53.03996 


-27.79847 


3.337 


P9 


3 


«i$s 


17.3/31 


<2.24E-17 


1.58E-16 


1.40E43 


26 


170 


53.04746 


-27.87047 


3.999 


P8 


4 


35±* 


308/359 


1.63E-15 


7.49E-15 


9.43E44 


1472 


197 


53.05771 


-27.93344 


6.071 


P8 


5 


Sfi+ 31 


88.2/124 


6.24E-16 


3.27E-15 


1.02E45 


134 


235 


53.07029 


-27.84564 


3.712 


S7 


1 


< 36 


33.0/27 


4.29E-17 


7.95E-17 


8.06E42 


27 


262 


53.07854 


-27.85992 


3.66 


SI 


1 


oc+22 
S£) -20 


123.6/165 


1.20E-16 


1.53E-15 


1.63E44 


208 


283 


53.08467 


-27.70811 


3.204 


P8 


1 


dJ -22 


58.0/80 


1.04E-16 


9.54E-16 


7.93E43 


77 


285 


53.08558 


-27.85822 


4.253 


P8 


3 


< 6 


29.7/31 


5.71E-17 


<9.13E-17 


1.24E43 


18 


331 


53.10271 


-27.86061 


3.78 


P10 


4 


< 23 


48.1/36 


4.61E-17 


<7.39E-17 


7.71E42 


22 


371 


53.11158 


-27.76789 


3.101 


P8 


4 


JO -10 


98.1/118 


1.22E-16 


7.70E-16 


5.70E43 


116 


374 


53.11204 


-27.86072 


3.724 


P9 


4 


01 n+284 
zid -114 


32.7/34 


<3.59E-18 


2.95E-16 


4.14E43 


22 


386 


53.11796 


-27.73439 


3.256 


S6 


1 


< 7 


62.6/64 


8.85E-17 


<1.42E-16 


1.06E43 


40 


403 


53.12196 


-27.93883 


4.762 


S2 


1 




66.9/73 


<5.46E-17 


1.30E-15 


2.68E44 


36 


412 


53.12442 


-27.85169 


3.7 


SI 


1 




208.0/225 


2.08E-16 


2.46E-15 


2.85E44 


373 


458 


53.13854 


-27.82128 


3.474 


S7 


1 


qo+74 


14.7/25 


1.59E-17 


2.52E-16 


2.64E43 


27 


485 


53.14658 


-27.87103 


7.620 


P8 


4 


406±^ 2 


74.2/94 


5.75E-17 


8.86E-16 


5.00E44 


87 


521 


53.1585 


-27.73372 


3.417 


P8 


4 


< 23 


47.6/50 


8.91E-17 


<2.40E-16 


2.05E43 


30 


528 


53.16158 


-27.85606 


3.951 


S5 


1 




81.6/111 


9.56E-17 


8.93E-16 


1.16E44 


126 


546 


53.16533 


-27.81419 


3.064 


SI 


1 


5211 


331.1/355 


6.75E-16 


6.26E-15 


4.73E44 


1072 


556 


53.17012 


-27.92975 


3.528 


P8 


3 


n 7 +9 
M '-10 


285.6/293 


7.59E-16 


1.07E-14 


1.14E45 


718 


563 


53.17442 


-27.86742 


3.61 


SI 


1 


g+2 
"-2 


257.6/307 


2.07E-15 


4.66E-15 


4.45E44 


1084 


573 


53.1785 


-27.78411 


3.193 


SI 


2 


q+2 


209.2/218 


8.10E-16 


1.65E-15 


1.19E44 


609 


588 


53.18467 


-27.88103 


3.471 


SI 


1 


< 2 


115.5/153 


6.45E-16 


1.03E-15 


8.41E43 


281 


642 


53.20821 


-27.74994 


3.769 


P9 


3 


< 13 


55.7/65 


9.17E-17 


<1.48E-16 


1.54E43 


18 


645 


53.20933 


-27.88119 


3.47 


S3 


2 


1511 


371.4/384 


3.16E-15 


1.03E-14 


9.24E44 


2053 


651 


53.21529 


-27.87033 


4.658 


P8 


3 


151±!1 


115.8/132 


1.33E-16 


2.38E-15 


4.64E44 


116 


674 


53.24004 


-27.76361 


3.082 


S3 


1 


< 7 


120.3/146 


6.95E-16 


1.39E-15 


8.56E43 


199 


700 


53.2625 


-27.86308 


4.253 


P8 


5 


"US 


73.8/90 


7.43E-16 


2.23E-15 


3.11E44 


102 


717 


53.28 


-27.79892 


4.635 


S4 


2 


87l£ 


47.2/73 


1.21E-16 


1.02E-15 


1.85E44 


41 


730 


53.29587 


-27.79317 


3.724 


S4 


1 




59.6/83 


<1.81E-16 


<1.46E-15 


1.65E44 


37 



(1) Identification number from lXue et alj j201lT ); (2) Right ascension and (3) declination from lXue et al.l j201l|); (4) adopted reds hift 
(see §|2]l; (5) referenc e for the adopted reds h ift - SI, S2, S3, S4, S5, S6 and S7: spectroscop i c redshift from|Szokolv et al. (2001), 
IVanzella et al. | 1120081 ). iPopesso ct al. (2009). Silverman et al. | l l2010l) , IVanzella et all ll201of),|Wuvts et al.l J2009l'l and E. Vanzella 



(private commun ication) , respectively; P8, P9 and P10: photometric redshift from lLuo et al. I ^OIOM Cardamonc ct al and 
Dahl en et alj d2010l . weighted solution), respectively; (6) criterion satisfied for source inclusion in the high-redshift sample (see §(2j: (7) 
best-fitting Nh as derived in S 13.21 in units of 10 22 cm -2 . We define Nh to be constrained if the lower limit (at the 90 per cent c.l.) on 
the best-fitting value is larger than zero, otherwise we report its upper limit. (8) best-fit value of Cash statistic over degrees of freedom; 
(9) soft-band (0.5 — 2keV) and ( 10) hard-band (2 — lOkeV) flux or upper limits from the best-fitting spectral model, in units of 
erg cm - 2 s — 1 ; (11) intrinsic (i.e. absorption-corrected) rest-frame 2 — lOkeV luminosity, in units of ergs -1 (typical errors range from 
~ 10 per cent for the brightest sources to a factor of ~ 2 for the faintest ones, fixing the redshift to the value reported in column 4); 
(12) full-band (0.5 — 7keV) net (i.e. background-subtracted) counts, as obtained from the spectral extraction. The full-band net-count 
number is only representative, being the energy range considered during the spectral fitting narrower, in most cases, in order to 

maximise signal-to- noise ratio (SNR). 
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Figure 5. Best-fitting Nh as a function of redshift for all the 
sources in the sample. Filled circles represent constrained values, 
while upper limits are reported as downward-pointing arrows. The 
dashed line connects the best-fitting Nh of XID 485 at z = 7.62, 
as reported in Tab. [Tl with the best-fitting value which would be 
derived if it were placed at z = 4.42 (open circle). See text for 
details. 



lLa Franca et all 120051 . lHasingerl |200St iBrusa et all l201(t> . 
possibly linked to radiative and mechanical feedback, sweep- 
ing away the obscuring gas and dust. We do not obtain a 
significant correlation between these two observables (see 
Fig- SI)- However, the luminosity range (43 ^ logZ/2-iokeV ^ 
45) is likely too narrow and the sample size too small to find 
such relation. 

Fig. [5] reports the best-fitting Nh plotted against the 
adopted redshift for each source. We note that the ma- 
jority (~ 68 per cent) of the sources are highly obscured 
(Nh > 10 23 cm~ 2 ). Furthermore, there is a hint for an in- 
creasing trend of column density with r edshift, as found in 
other works (e.g. iTreister fc Urrv||2006l ). However, a possi- 
ble bias due to the high redshift and low counting statistics 
can be present and will be investigated in the next section. 
The trend could be driven by the three sources at z > 5, 
which have the most uncertain redshift information (see 
§ [2j) . If they were placed at z = 3.7 (the median redshift 
of the sample), the best-fitting Nh would be 50t.ll x W 22 , 
42±H x 10 22 and 118±^ x 10 22 cm -2 for XID 139, 197 and 
485, respectively (i.e., they would still be in the highly ob- 
scure d class). Assuming t he photometric redshift (z = 4.42) 
from ISantini et all (|2009T l for XID 485, the only source m 
the sample at z > 5 for which more than one photometric 
redshift is provided by the catalogues that we considered, 
the best-fitting column density would be I45I3 2 x 10 22 cm -2 
(Fig. [5l open circle). We recall that the low photon counting 
statistics allow us to assume only a simple spectral model, 
which doe s not take into acco unt complex absorption struc- 
tures (e.g. lBianchi et al"] |2012). which could affect the results 
on the column density. 

The four sources with Nh > 10 24 cm~ 2 were also fitted 



with the XSPEC model plcabs (power-law emission trans- 
mitted through a spherical distribution of cold an d dense 
matt er, taking into account Compton scattering; lYaqoobl 
Il997f) . which reproduces X-ray emission from Compton- 
thick sources better than a simple absorbed power-law, but 
no significant change in the best-fitting spectral parameters 
was found. 

We also investigated the possible presence of a neutral 
iron Kq line in the spectra by adding a narrow (a = 0.01 
keV) Gaussian component (zgauss model in XSPEC) at 
6.4 keV rest-frame. We could only find upper limits to the 
rest-frame equivalent width (EW) for most of the sources, 
of the order of a few 100 eV on average. This result is 
easily explained given that the observed equivalent width 
scales with (1 + z)^ 1 and therefore only loose constraints 
can be placed to the line in spectra of high-redshift AGN 
(especially if the signal-to-noise ratio and the number of 
net counts are low and the accuracy in the redshift mea- 
surement is limited). The only source for which the iron 
line was easily visible in the spectrum a nd statistically 
significant is XID 412 (i.e. CDF-S 202 in iNorman et all 
2002). In this case, we allowed the rest-frame energy of 
the line to vary during the spectral fitting in a relatively 
large range. The best-fitting rest-frame energy, EW and Nh 
that we derived for this source (E = 6.53 keV, EW = 
738t 4 ^ eV, N H = 91±* 4 x 10 22 cm" 2 ) are in good agree- 
ment with those found by p revious Chandra (INorman et all 
120021 ) and XMM-Newton l|Comastri et all l201ll ) observa- 
tions. The soft (Fb. 5 _ 2 keV = 1-76 x 10" 16 erg cm" 2 s _1 ) and 
hard (F2-iokeV = 2.41 x 10~ 1 5 ergcm~ 2 s~ 1 ) flu xes are con- 
sistent with those reported bv lTozzi et alj (120061 ) using the 1 
Ms CDF-S data, although they assumed a different spectral 
model (r = 1.8 and Nh = 1.5 x 10 24 cm~ 2 ). 



3.3 Nh distribution 

The distribution of the best- fitting Nh (Fig. O derived from 
spectral analysis, assuming an absorbed power-law with 
fixed F = 1.8 as fitting spectral model, strongly suggests 
the presence of an important fraction (ps 68%) of highly ob- 
scured sources (Nh > 10 23 cm -2 ). This result could be af- 
fected by a systematic Nh overestimate due to the high red- 
shift and low counting statistics, characteristic of the sam- 
ple: the photoelectric cut-off in spectra of mildly-obscured 
(10 22 < N H < 10 23 cm" 2 ) sources is redshifted out of the 
Chandra band-pass already at z = 3 — 4, making very diffi- 
cult the search for the best-fitting Nh- Indeed, this type of 
spectrum might be indistinguishable from a spectrum of an 
unabsorbed source at high redshift. Furthermore, statistical 
fluctuations, enhanced by the low fluxes, are expected in the 
spectra and could simulate a spurious photo electric cut-off 
at high rest-frame energies (|Tozzi et alj|2006h . 

In order to estimate the probability (Pij) that the X-ray 
spectrum of a source with intrinsic Nh in a given range j 
is fitted with a best-fitting column density constrained in a 
range i, we made extensive use of simulations. As in Tab.[TJ 
we define Nh to be constrained if the 90% c.l. lower limit 
on the best-fitting value is larger than zero. We divided the 
sample in 5 Nh bins, as reported in Tab. [5] The number of 
sources Ni observed in each absorption bin i is then given 
by: 
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Figure 6. Observed distribution of the best-fitting Nh- The 
black line represents the constrained values; the red line with 
leftward pointing arrows represents the upper limits. 



N, 



(1) 



where Xj is the number of sources with intrinsic Nh in bin 
j and Ni is the number of sources observed with best-fitting 
Nh in bin i, for i, j = A, B, C, D, E. In this framework, Ni 
is known (black histogram in Fig. [6} and Pij can be derived 
from simulations. 

We simulated 1000 spectra with T — 1.8 and 100 net 
counts (close to the median value of our sample), intrinsi- 
cally obscured by a column density in the j-th Nh bin (for 
each bin we used the value reported in the third column of 
Tab. [2|, considering sources at z — 4 (similar to the char- 
acteristic redshift of our sample), and simply counted how 
many times the simulated spectra are fitted with a best- 
fitting Nh constrained in the i-th bin (see Appendix A for 
the detailed description of the simulation procedure). Bins 
A, B and C are indistinguishable at z > 3 and hence are 
merged into a single bin (representing unobscured or mildly- 
obscured sources with Nh < 10 23 cm -2 ), hereafter referred 
to as (ABC). 

The value of each Pij is shown in Tab. [3] We immedi- 
ately note that no spurious effect seems to be important in 
bins D and E (i.e. spectra of intrinsically strongly obscured 
sources are efficiently fitted by a model obscured by a col- 
umn density within the correct Nh bin) while, as expected, 
the best-fitting Nh is overestimated in most cases (ss 64%) 
when the intrinsic value is in the bin (ABC). 

By inverting eq. Q), the intrinsic Nh distribution can 
now be derived (Fig. [7|: :E(abc) ~ 5.5, xo ~ 15.5 and 
xe = 4, to be compared with the observed distribution: 
N(abc) — 2, No = 19 and Nb = 4. If we conservatively 
count all the 9 sources for which only upper limits could be 
obtained (red line in Fig.[5| in the first bin, the resulting Nh 
distribution is the one shown in Fig. [8] W e adde d the pre- 
dictions of the iGilli. Comastri fc Hasingeij (|2007l . hereafter 



Table 2. 


Nh bins used in 


simulations 


i,j range: log(^^.) 


sim: log( Nh . > ) 

ov cm z ' 


A 


< 21 


0.0 


B 


21-22 


21.5 


C 


22-23 


22.5 


D 


23-24 


23.5 


E 


>24 


24.5 



Sources in the i-th or j-th Nh bin (first column), which accounts 
for the Nh range shown in the second column, are simulated 
assuming the column density reported in the third column. 



Table 3. Final probability factors Pij that a source with intrinsic 
Nh in a given bin j is observed in a bin i. These are derived from 



simulations 


using r = 1.8 


and T 


= 1.6. 


p 


i=(ABC) 


i=D 


J'=E 


r = i.8 


i=(ABC) 


0.361 


0.0 


0.0 


i=D 


0.639 


1.0 


0.0 


i=E 


0.0 


0.0 


1.0 


r = i.6 


j=(ABC) 


0.147 


0.0 


0.0 


i=D 


0.853 


1.0 


0.0 


i=E 


0.0 


0.0 


1.0 



GCH07) X-ray background (XRB) synthesis model com- 
puted with the pompaIII tool for each Nh bin, considering 
sources with luminosity 42 < Log °- 5 - 2 _ k 1 ° v <j 47 in a red- 
shift range between z m in = 3 and z max — 8. A high-redshift 
exponential decline in the space density of AGN, similar to 
that of bright quasars l|Gilli et al.ll2011f) . is also considered. 
The assumption of the decline affects only the normalization 
and not the shape of the distribution. The predictions are 
corrected for the sky coverage (i.e. the sky area in physi- 
cal units covered by the survey at different fluxes, see next 
section), to take into account the smaller sky area covered 
by the survey at low fluxes. Obscured sources have a flux 
lower than unobscured or less obscured ones with the same 
luminosity at the same redshift and, hence, can be detected 
over a smaller, deeper area. 

The Nh distribution, corrected for the spurious overes- 
timate of the column density, shows a number of strongly 
obscured sources that is slightly larger than the predictions 
of the model we considered. Conversely, the number of un- 
obscured or mildly-obscured sources is smaller than the pre- 
dicted number. However, we find an agreement between the 
observations and the predictions at the 2a confidence level 



7 POrtable Multi-Purpose Application for AGN counts 
| |http://www. bo. astro.it/~g illi / counts . html ) 
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Figure 7. Intrinsic Nh distribution (solid line) with statistical 
errors, derived correcting the observed one (dotted line, plotted 
in the same bins, with statistical errors slightly shifted for visual 
purpose), resulting from spectral analysis assuming V = 1.8, with 
the probability factors obtained from simulations (see § 13,31 and 
Appendix |A]| . Only constrained values of Nh are considered. 



(the x 2 test returns a value of 6.0 for 3 degrees of freedom 
and no free parameter). 

Following the results of § 13.21 where, considering an ab- 
sorbed power-law with F free to vary as starting spectral 
model, we derived a mean photon index < V >~ 1.6, we in- 
vestigated the effects of assuming a flatter photon index on 
the column density distribution. Therefore, we fitted again 
the X-ray spectra of all sources in the sample using the same 
spectral model with photon index fixed to T = 1.6. The re- 
sulting Nh distribution was corrected for spurious overesti- 
mates using the probability factors (Tab. [3| obtained from 
simulations, similarly to the T = 1.8 case, but now using 
r = 1.6 to simulate and fit the spectra. Red points in Fig. [8] 
represent the intrinsic Nh distribution in the F = 1.6 case. 
We derived upper limits on Nh for two sources (XID 132 and 
700, which were previously counted in bin D), while their 
best-fitting column density was constrained using V = 1.8. 
This result produces a different number of sources in the 
two shaded histograms in Fig. [8] Therefore, the number of 
objects in each bin of the corrected Nh distribution is, in 
this case, £(abc) ~ 13.6, id ~ 5.4 and ie = 4. If we place 
all the 11 sources with an upper limit on Nh in the first 
bin, as in Fig. [8] the distribution is fully consistent within 
2(7 with the prediction of the model (\ 2 = 4.75 with 3 de- 
grees of freedom and no free parameter). We underline that 
the model predictions are computed with the same level of 
incompleteness as the observed sample. 
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Figure 8. Intrinsic Nh distribution assuming T = 1.8 (black 
points) and T = 1.6 (red points, which are slightly shifted from 
the centre of the bins for visual purposes). The first bin includes 
sources with Nh < 10 23 cm _ 2 (i.e. consistent with no or "moder- 
ate" absorption) and all the sources for which only upper limits 
on Nh were obtained (9 and 11 sources in the two cases, repre- 
sented as grey and orange shaded areas, respectively). The third 
bin contains sources with Nh > 10 24 cm — 2 . Poissonian errors are 
computed for the number of sources in each bin. The histogram 
is the predictions computed with the POMPA tool and corrected 
for the sky coverage of the survey (see text). The predictions are 
not normalized to the observed number of sources. 



3.4 The logN-logS of the high-redshift sample 

We derived the logN-logS of the high-redshift sample in the 
soft band, where the survey sensitivity is the highest, folding 
the number TV of sources with flux larger than S with the 
survey sky-coverage at that flux. We assumed the soft fluxes 
reported in Tab. [T] derived from our spectral analysis. The 
sky-coverage was obtained by multiplying the nominal sky 
area covered by the survey (« 0.13 deg 2 ; IXue et al J 120111 ) 
and the fraction of the f.o.v. covered at different fluxes, de- 
rived ap plying the CIAO tool DMIMGHIST on the sensitiv- 
ity map l|Xue et al.ll201ll ) , which assigns to each point of the 
f.o.v. the flux limit of the survey at that position. In order 
to be consistent with the computation of the sky coverage 
(derived in the 0.5 — 2keV band), we excluded four (w 12% 
of the total sample) sources which were not detected in the 
soft band (see Tab.rjJ. 

In Fig. [9] we c omp are our data w it h tho se provided by 
iFiore et all (f2012h and iLehmer et all (l2012h for the 4 Ms 
CDF-S and with th e XMM l|Brusa et alj|2009l ) and Chandra 
dCivano et aLlbOllT ) COSMOS data. The predictions of the 
IGCH07 XRB synthe sis model (solid and dashed lines) and 
the lAird et~al] feoich X-ray luminosity function (X LF) evo- 
lutionary model (dotted line) are also shown. The GCH07 
model is based on a luminosity dependent density evolu- 
tion (LDDE) model of the XLF with redshift, and it al- 
lows the presence of a high-redshift (z > 2.7) de cline in the 
AGN space density to be taken into account. The lAird et al.l 
(2010) model assumes a luminosity and density evolution 
(LADE) of the XLF. 

Our results for the z > 3 sam ple (top panel) are in 
good agreement with the IGCH07I model predictions as- 
suming an exponential decline in high-redshift AGN space 
density (continuous line), down to flux .Fo.5-2keV ~ 4 x 
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Figure 9. Binned LogN-logS of the high-redshift sample, consid- 
ering all sources (z > 3, top panel) and only those at z > 4 (bot- 
tom panel), with associated Poissonian errors . Black filled circle s 
are from this work, green triangles are from iFiore et al 
blue squares are fro m lCivano et alj j201ll ) , red pent agons are from 
iBrusa et al. and purple hexagons are from lLehmer et al.l 

l|2012l l. Empty black circles in the bottom panel represent our 
data if the three sources at the highest reds hifts were placed at 
3 < z < 4 (see text). The predictions of the IGCH07I XRB syn- 
thesis model, computed with the POMPA tool, either assuming 
(solid line) or not (dashed line) an ex ponential decline in the 
AGN space density at high redshift, and lAird et alj i201Gt ) model 
are also shown. The yellow shaded area represents the regions in 
which the LogN-LogS would be placed if more conservative or 
shallower selection criteria were used as well as the grey area, 
which also includes the 39 sources with no redshift in the consid- 
ered catalogues (see §dj. 



10 — 17 ergcm -2 s -1 , a factor « 10 fainter with respect to pre- 
vious surveys (i.e. down to Lx ~ 10 43 ergs -1 , the minimum 
luminosity detectable at z > 3 in the 4 Ms CDF-S): we ex- 
tend to lower fluxes the re sults obtained in t h e CO SMOS 
field by IBrusa et~ai1 (|2009l ) and ICivano et all (|201ll ), con- 
firming an evolution at high redshift similar to that of the 
more luminous AGN |Gilli et al.ll201lh . 

Considering only the z > 4 subsample (Fig. [9] bot- 
tom panel, filled circles), our data lie above the model 
predictions. We include here the presence of 3 sources at 
5 < z < 7.7, whose redshifts are determined on the ba- 
sis of relatively uncertain photometric information. If these 
sources were placed at 3 < z < 4, a good agreement would 
be obtained (Fig. [9l bottom panel, open circles). 

The scenario without the decline (Fig. [5] dashed lines) 
is excluded by our data both at z > 3 and z > 4. 



4 DISCUSSION 

The Nh distribution derived from spectral analysis, cor- 
rected for the overestimate of the column density, suggests 
the presence of a significant fraction (~ 57 per cent) of 
strongly obscured (N H > 10 23 cm~ 2 ) AGN at z > 3. This 
result was obtained by fitting the X-ray spectra with an 
absorbed power-law with F = 1.8, i.e. the t ypical photon 
index observed in larg e AGN samples (e.g. iTurner et al.l 
1 19971 ; iTozzi et al J 12009 ). Adopting F = 1.6, we derived a 
slightly different intrinsic Nh distribution, with more un- 
obscured sources, with respect to the previous case. Since 
the current data are sampling a rest-frame energy range up 
to at least Ei, lt ~ 20 keV, a hardening of the X-ray spec- 
tra due to reflection might be expected at observed ener- 
gies E bs ~ 4 — 5keV. However, the low photon-counting 
statistics prevented us from using relatively complex spec- 
tral models and the T = 1.6 index may simply incorporate 
the hardening due to reflection. Therefore, the "true" Nh 
distribution is expected to lie between the two distributions 
derived assuming V = 1.8 and V = 1.6 (Fig. |HJ). To further 
investigate the possible impact of a reflection component, we 
computed the stacke d, rest frame, X-ray spectra (using the 
method described in llwasawa et al]|2012l ) dividing the sam- 
pie in four N H bins ( 10 / c *_. < 1.3; 1.3 < 1f , a3 ^_ a < 6; 

N 



6 < 1<s *?? m -* < 9; 10 ^_ 2 > 9; see table 1). The result- 
ing best fit slope, obtained by a joint fit of the the various 
Nh bins and representing the intrinsic continuum above 10 
keV, is r ~ 1.69^q'q7- This slope is suggestive of a lim- 
ited contribution from reflection and is well within the val- 
ues adopted to compute the absorption distribution (Fig. 
8). Independently from the assumed spectrum, the data 
are in agreem ent at th e 2a confidence level with the pre- 
diction of the IGCH07I XRB synthesis model, assuming an 
exponential decline in the AGN space density at z > 2.7 
(Fig. |HJ. The best-fitting column densities of the 19 high- 
luminosity (Z/2-iokcV > 10 44 ergs -1 ) sources are well con- 
strained at high values (~ 84 per cent of them has nominal 
Nh > 10 23 cm -2 ). Applying the correction factors derived 
assuming F — 1.8 (see § 13. 3[) . we derived that the intrin- 
sic fraction of obscured sources (A^sr > 10 23 cm~ 2 ) in this 
subsamp le is ~ 66 per cent. 

The iGCHOTl model assumes that the fraction of ob- 
scured AGN does not evolve with redshift. A possible evolu- 
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tionary scenario emerging from recent works suggests that 
this behaviour could be valid o nly for qua s ars (Q SOs) with 
^2-10 kcV sC lO^ergs" 1 (see iGilli et all |2010| ). which is 
the median luminosity of our sample. Indeed, evidences 
are found that the activity of luminous (Z/2-iokcV ^ 
10 44 ergs _1 ) AGN, triggered by major gas-rich mergers, 
occurs in short bursts (with timescal es of the order of 
~ 0.01 Gyr; e.g. I Alexander et al. | 2005). In this case, the 
gas accretion is chaotic ( Hopkins et alJ 20081 ) and may pro- 
duce a large covering factor and column densities. Since 
the merging rate increases with redshift, the fraction of ob- 
scured AGN may then increase as well , as reported in e.g. 
lLa Franca et al.l |2005h ■ Frlasingerl (|2008h and Iwasawa et al. 
(2012, submitted). At lower luminosities, the accretion onto 
the SMBH is probably driven by secular processes (~ 1 Gyr; 
e.g. iDaddi et al.ll200?1 ; lElbaz et al.ll201ll ) and the accretion 
is expected to be smoother and symmetrical. Therefore, the 
fraction of obscured AGN would be quite constant with red- 
shift, given that obscuration would be a function of the sys- 
tem geometry only. Our results fit well in this framework. 

Fig. [5] shows the number counts derived in this work 
compared to previous determinations from the literature 
and theoretical models. The CDF-S source counts down to 
-Fo.5-2k c v ~ 4 x 10 _17 ergcm~ 2 s~ 1 lie on the XRB model if 
an exponential decline of the A GN X-ray luminos ity function 
is as sumed, similarly to XMM (|Brusa et al . 2009) and Chan- 
dra (|Civano et al.ll201ll) COSMOS results at brighter fluxes. 
In particular, we note the agreement at bright fluxes with 
the results from COSMOS surveys, for which a major ef- 
fort w as made to obtain spectroscopic and photometric red- 
shifts jSalvato et al.ll2009l ; iBrusa et all [20091 ; ISalvato et all 
l201ll . and re ferences therein). Our results are consistent 
also with the lAird et al.l (|2010l ) model if the least conser- 
vative bounds (grey area in Fig. [9| are assumed (see be- 
low for the computation of the error budget). At fluxes 
2 x 10~ 16 sC F . 5 _2kov «S 2 x lO^ergcm-V 1 (i.e. for 
the brightest sources in the high-redshift s ample) t he pre- 
diction of this model is very similar to the GCH07 model, 
but at fainter fluxes the t wo mode ls deviate significantly. 
Our results better trace the GCH07 m odel, w hile data from 
iFiore et al. I ||2012|) an d iLehmer et all (|2012l ) are placed on 
the lAird et all (|20ld ) model. 

There is consiste ncy in the numb er co unts that we 
presen t with data from lFiore et al.l |2012T ) and lLehmer et all 
l|2012l ) down to fo.5-2kcV~10 erg cm s , while our re- 
sults are a factor ~ 2 — 3 lower at fainter fluxes. This is not 
surprising considering that different selection techniques are 
applied, ours being more conservative than the ones used in 
those works. Indeed, the logN-logS in lFiore et al.l (|2012T ) is 
derived for 54 optically selected sources at 2 > 3, using spec- 
troscopic and photometric redshifts and drop-out techniques 
and applying a new detection algorithm on the 4 Ms CDF-S 
observations, which made possible to reach a flux limit of 
Jb.5-2k eV ~ 1-2 x 10~ 17 erg cm~ 2 s~ 1 . The samples consid- 



ered in lLehmer et al.l |2012l ) and in this work are collected 
from the same cat a logue of X-ray selected sources. How- 
ever. ILehmer et al.l (|2012l ) assumed the redshift adopted in 
IXue et al l (|201ll ). while we collected additional information 
from other catalogues and used a stricter selection proce- 
dure (as described in §H]). There fore, we rejected a nu mber 
of sources which are included in ILehmer et al.l {2012), and 
we were able to collect redshift information for some objects 



which have no redshift in IXue et all (|201lh . as in the case 
of XID 331. We wou ld obtain results largely consistent with 
ILehmer et all (|2012T ) using their sample. 

To compute the "error budget" on the number counts, 
we investigate how many sources would be included in 
the sample if we relax the third, fourth and fifth crite- 
ria described in § [2] Relaxing the third criterion and in- 
cluding in the final sample soft X-ray sources for which 
N p hot,>3 — N p hot,<3 ^ 0, the net gain in the final sample 
would be of 4 sources. As for the fourth criterion, 14 sources 
were excluded from the sample because only one of the two 
available photometric redshift is at 2 > 3. We consider these 
objects to have a 50% chance of being at z > 3, and ran- 
domly include seven of them in the error budget. Finally, 
if we relax the fifth criterion and include in the sample all 
sources with only one photometric redshift, with the con- 
dition of being at z > 3, regardless of its error, we would 
add 4 sources to the final sample. If all these changes were 
adopted, the sample would consist of 49 sources. 

The gold shaded area in Fig. [9] indicates the regions in 
which the logN-logS would lie considering the error bud- 
get 0. For the z > 3 number counts (upper panel), the 
lower bound is obtained if only sources selected via the first, 
second and third criteria of § [2] are included in the high- 
redshift sample (i.e. sources with spectroscopic redshift or 
N p hot,>3 — N p hot,<3 ^ 3). This is clearly a very conservative 
assumption. The upper bound is obtained if we relax the 
third, fourth and fifth criteria, as discussed above, adding 
15 sources in total. 

Four out of these 15 sources would be at z > 4: two 
would be selected relaxing the fifth criterion and two relax- 
ing the fourth criterion. As previously described, we ran- 
domly assumed half of the latter sources to be at z > 4. The 
resulting 3 sources were added to the z > 4 sample to com- 
pute the upper-bound of the gold shaded area (bottom panel 
in Fig. [5]), which accounts for the error budget of the z > 4 
number counts. The lower bound is obtained, as before, con- 
sidering only the first, second and third criteria, and consists 
of one source with a spectroscopic redshift z > 4 (XID 717, 
while we discarded XID 403 because it is not detected in 
the soft band) and 2 sources selected through the third cri- 
terion. As expected, because of the low number of sources 
detectable at z > 4, only loose constraints can be placed on 
the LogN-LogS. Larger samples of X-ray selected AGN at 
2 > 4 would be needed to draw more solid conclusions. 

Thirty-nine sources detected in the soft band in 
IXue et al.1 l |201ll ) have no redshift in any catalogue that we 
considered. A not negligible fraction of them are probably 
AGN at z > 3. The grey shaded regions in Fig. [9] show 
how the gold area would enlarge if all these sources would 
be accounted for. Therefore, these areas represent the hypo- 
thetical and unlikely cases in which all the 39 sources were 
at z > 3 (2 > 4). 



8 We assigned the soft X-ray fluxes reported in lXue et alJ ll201l[ ) 
to the sources not included in the high-redshift sample for the 
computation of all shaded areas. 
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5 CONCLUSIONS 

The detection of a significant number of z > 3 AGN down 
to Lx ~ 10 43 ergs -1 is made possible by t he faint flux limi t 
reached in the 4 Ms CDF-S observations jXue et al.ll201lh . 
Moreover, X-rays are less biased against obscuration with 
respect to the optical band. The large multi-wavelength cov- 
erage of the CDF-S allows us to retrieve spectroscopic and 
photometric redshifts for most of the sources in the field. 
In this framework, physical and evolutionary properties of a 
relatively large sample of high-redshift AGN can be studied. 

The sample consists of 34 sources in the range 3 < z ^ 
7.6, with a median redshift of z — 3.7. About 45 per cent of 
them (15 sources) have a spectroscopic redshift. 

X-ray spectra of all the sources were extracted from 
the 4 Ms CDF-S observations and fitted with an absorbed 
power-law spectral model with photon index fixed to T = 
1.8. The low photon-counting statistics prevents us from as- 
suming a more complex spectral model: the median full- 
band (0.5 — 7keV) net-count number is 80. The median 
rest-frame absorption-corrected luminosity of the sample is 
-I/2-io kov « 1.5 x 10 44 ergs _1 . 

We summarise the conclusions reached in this work as 
follows: 

(i) The observed column density distribution is strongly 
peaked between 23 < log( N 1 2 ) < 24. The low number 
of counts and the high redshift of the source sample may 
lead to overestimate the intrinsic absorption. To account for 
this effect, we ran extensive X-ray spectral simulations and 
derived correction factors to apply to the observed distribu- 
tion. Therefore we obtained an estimate of the intrinsic Nh 
distribution. The results are consistent with a non-evolving 
obscured fraction of AGN with L x w 10 ergs -1 with red- 
shift within 2cr. 

(ii) The number counts of our z > 3 AGN sample are 
consistent with a decline in the AGN space density at 
high redshift, extending at faint soft fluxes (-Fo.5-2k e v ~ 
4 x 10~ 17 ergcm ~ 2 s~ 1 ) the behavio ur de termined at b r ighter 
fluxes found bv lBrusa etail (120091 ) and lCivano et all (|201ll ) 
in the COSMO S field. We u s ed a m ore conservative approac h 
with respect to lFiore et ail (|2012f ) and lLehmer et all <|2012h : 
if we significantly relax our selection criteria, our results are 
in a relatively good agreement at bright fluxes with their 
data, while at the faint end the discrepancy is a factor ~ 2. 

These results would be improved by increasing the size 
and the reliability of the sample (e.g. gathering additional 
reliable redshifts via ultradeep optical-IR spectroscopic cam- 
paigns) and the quality of the X-ray spectra (e.g. increas- 
ing the Chandra exposure in this field). The accuracy of 
the photometric redshifts would be significantly improved 
by adding the deep Cosmic Assembly N ear-IR Deep Legacy 
Survey (CANDELS) H and J-band data (|Grogin et alj|201ll : 
iKoekemoer et al.ll201ll ). 
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APPENDIX A: DERIVING THE CORRECTION 
FACTORS FOR THE OBSERVED N H 
DISTRIBUTION FROM SIMULATIONS OF 
X-RAY SPECTRA 

The correction factors discussed in § 13.31 are derived from 
spectral-fitting simulations using XSPEC, assuming five dif- 
ferent degrees of obscuration (see Ej |3,3l and Tab. [2}. The five 
cases are simulated by runs of 1000 simulations each. The 
spectra are simulated at z — 4 with an input normalization 
tuned to obtain a distribution of net counts peaked at « 100 
(with a resulting dispersion of a w 11%), using the response 
files of a real source at an off-axis angle 6 — 4arcmin. The 
results do not vary significantly if different response files are 
used, but they are sensitive to the number of net counts and 
redshift. To be consistent with the real cases, for these pa- 
rameters we assumed values close to the average ones of the 
sample. 

The simulation procedure can be summarized in the 
following steps: 

(i) Definition of an input model: we considered a power- 
law with photon index fixed to T — 1.8 as input model, ab- 
sorbed by five different column densities, whose values vary 
between the five considered cases (Tab. [2] third column). 
Galactic absorption (Nh — 7 x 10 19 cm -2 ) is also included. 
The desired number of counts is obtained by properly ad- 
justing the normalization of the power-law. 

(ii) Simulation of the spectrum: the fakeit command in 
XSPEC can simulate a spectrum with background, given a 
starting model. It requires as input the response and an- 
cillary files, a real background file (with the POISSERR 
keyword set to TRUE) and the desired exposure time. We 
set it to 3.5 x 10 6 s to be consistent with the average real 
case. The output consists of a fake source and a fake back- 
ground spectrum file. The same BACKSCAL keyword, read 
from the header of the input background file, is assigned 
by the fakeit procedure to both simulated spectra, which 
are therefore assumed to be extracted from regions with 
an equal area. This approach is not correct and causes the 
source and background counts to be wrongly scaled (i.e. a 
strong underestimate of the source net-count rate). To ad- 
dress this problem, before simulating the spectra we multi- 
plied the EXPOSURE keyword of the real background file 
by a factor equal to the ratio between the BACKSCAL pa- 
rameters of the real background and associated source file. 
We checked that the errors on the net-count rate were cor- 
rectly computed (i.e. propagating the errors of the total and 
background-only count-rates) . 

(iii) Spectral fitting: after having grouped the simulated 
spectra, using the FTOOL grppha, with at least 1 count 
per bin, and assigned to them the response files and the 
fake background spectra, the spectra are repeatedly fit in 
the energy range E = 0.5 — 5keV with an absorbed power- 
law model with Nh free to vary and photon index fixed to 
r = 1.8, accounting also for the Galactic absorption (fixed 
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to Nh = 7 x 10 19 cm~ 2 ). Cash statistics were employed in 
the fitting. The fitting step (which use a local minimization 
algorithm) is alternated to the computation of the errors of 
the various free parameters (which is a more global opera- 
tion) for a few times, in order to avoid the fit to be stuck in 
a local minimum. 

For each of the five considered cases of input Nh, we 
obtained as output 1000 best-fitting values of Nh- Following 
the procedure adopted to fit the real spectra, we define Nh 
to be constrained if the lower limit (at the 90 per cent c.l.) 
on the best-fitting value is larger than zero, otherwise we 
consider its upper limit (e.g. if the fit returns a best-fitting 
value N H = 10l^ x 10 22 cm~ 2 , we consider N H < 22 x 
10 22 cm" 2 ). 

The probability pij that a simulated source in the j- 
th Nh bin is observed in the i-th Nh bin, are derived by 
counting the number of times a constrained best-fitting Nh 
falls in bin i and normalized this number to 1000. In order to 
conserve the number of sources when the correction factors 
are applied on the real distribution, the ptj factors must 
be rescaled to the total probability to constrain Nh (i.e. 
we excluded all cases in which only upper limits could be 
found): 



P — 



Pij 



EE 



(Al) 



As discussed in § 13.31 bins A, B and C are merged into a 
single bin (ABC), since they are indistinguishable at z > 3. 
Therefore, the probability factors, -P;(abc)i are: 



(ABC)(ABC) 



for 



(A2) 



PiA+PiB+PiC 



and 



3 = A, B, C; P i( abc) - 
P(abc) . = P A]+ P B]+ Pc 3 for . f . = D and E 

The same procedure is then repeated using T — 1.6 
instead of F = 1.8. The Pij factors, for i and j = (ABC), D 
and E, are those used in § 13.31 



APPENDIX B: 
SAMPLE 



X-RAY SPECTRA OF THE z > 3 
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Figure Bl. X-ray spectra of the 34 sources in the high-redshift 
sample, fitted with an absorbed power-law model with V = 1.8 
(see § 13.211 . The ratio of data to folded model is shown and the 
spectrum is binned (at ^ 3<r detection signifi cance) for displa y 
purposes only. Identification numbers are from lXue et alj J201ll) . 
The spectroscopic or photometric redshift is reported for each 
source, as well as the best-fitting column density, accordingly to 
Tab.Q] (Continue) 
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Figure Bl. (Continue) 
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XID285 zphol=4.253 
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High-redshift AGN in the CDF-S 
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